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Abstract:  We  investigated  molecular  interactions  involved  in  the  selective  binding  of  several  short  peptides 
derived  from  phage-display  techniques  (8-12  amino  acids,  excluding  Cys)  to  surfaces  of  Au,  Pd,  and 
Pd-Au  bimetal.  The  quantitative  analysis  of  changes  in  energy  and  conformation  upon  adsorption  on  even 
{111}  and  {1 00}  surfaces  was  carried  out  by  molecular  dynamics  simulation  using  an  efficient  computational 
screening  technique,  including  1000  explicit  water  molecules  and  physically  meaningful  peptide  concentra¬ 
tions  at  pH  =  7.  Changes  in  chain  conformation  from  the  solution  to  the  adsorbed  state  over  the  course  of 
multiple  nanoseconds  suggest  that  the  peptides  preferably  interact  with  vacant  sites  of  the  face-centered 
cubic  lattice  above  the  metal  surface.  Residues  that  contribute  to  binding  are  in  direct  contact  with  the 
metal  surfaces,  and  less-binding  residues  are  separated  from  the  surface  by  one  or  two  water  layers.  The 
strength  of  adsorption  ranges  from  0  to  -100  kcal/(mol  peptide)  and  scales  with  the  surface  energy  of 
the  metal  (Pd  surfaces  are  more  attractive  than  Au  surfaces),  the  affinity  of  individual  residues  versus  the 
affinity  of  water,  and  conformation  aspects,  as  well  as  polarization  and  charge  transfer  at  the  metal  interface 
(only  qualitatively  considered  here).  A  hexagonal  spacing  of  ~1 .6  A  between  available  lattice  sites  on  the 
{1 1 1}  surfaces  accounts  for  the  characteristic  adsorption  of  aromatic  side  groups  and  various  other  residues 
(including  Tyr,  Phe,  Asp,  His,  Arg,  Asn,  Ser),  and  a  quadratic  spacing  of  ~2.8  A  between  available  lattice 
sites  on  the  {100}  surface  accounts  for  a  significantly  lower  affinity  to  all  peptides  in  favor  of  mobile  water 
molecules.  The  combination  of  these  factors  suggests  a  “soft  epitaxy”  mechanism  of  binding.  On  a  bimetallic 
Pd-Au  {111}  surface,  binding  patterns  are  similar,  and  the  polarity  of  the  bimetal  junction  can  modify  the 
binding  energy  by  ~10  kcal/mol.  The  results  are  semiquantitatively  supported  by  experimental  measurements 
of  the  affinity  of  peptides  and  small  molecules  to  metal  surfaces  as  well  as  results  from  quantum-mechanical 
calculations  on  small  peptide  and  surface  fragments.  Interfaces  were  modeled  using  the  consistent  valence 
force  field  extended  for  Lennard-Jones  parameters  for  fee  metals  which  accurately  reproduce  surface  and 
interface  energies  [Heinz,  H.;  Vaia,  R.  A.;  Farmer,  B.  L.;  Naik,  R.  R.  J.  Phys.  Chem.  C  2008,  112, 
17281-17290], 


1.  Introduction 

A  variety  of  inorganic— biological  nanostructures  and  bio- 
mimetic  materials  with  distinct  functionalities  can  be  assembled 
using  nanoparticles  of  various  chemical  composition,  peptide 
linkers,  surfactants,  and  polymer  matrices. 1-16  Peptides  may 
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regulate  crystal  growth  of  the  inorganic  material, 1,3-8,13  attract 
metal  nanoparticles  of  a  certain  composition  and  size  to  form 
characteristic  metal-coated  wires,3-5  assemble  in  predefined 
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Interactions  of  Peptides  with  Au,  Pd,  and  Pd-Au  Surfaces 
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Table  1.  Peptide  Names,  Binding  Targets,  and  Amino  Acid 
Sequences 


name 

binds  to 

peptide  sequence" 

A3 

Au,  Pd" 

+Ala-Tyr-Ser-Ser-Gly-Ala-Pro-Pro-Met-Pro-ProPhe_ 

Fig— Na3 

Au,  Pd" 

+Asp_-Tyr-Lys+-Asp_-Asp_-Asp_-Asp_-Lys+_  •  3  Na+ 

Pd2— Cl 

Pd* 

+Asn-Phe-Met-Ser-Leu-Pro-Arg+-Leu-Gly-His-Met-His-  • 

cr 

Pd4— Cl 

Pd* 

+Thr-Ser- Asn- Ala- V al-His-Pro-Thr-Leu- Arg+-His-Leu-  • 

cr 

Flgd— Na2 

+Asp--Tyr-Lys+-Asp--Asp--Asp--Lys+-  •  2  Na+ 

GlylO 

Glyio 

ProlO 

Pro10 

“  References  7,  8.  b  These  peptide  sequences  from  phage  display  are 
newly  reported  here.  c  The  approximation  of  the  charge  state  equivalent 
to  pH  =  7  in  the  models  is  indicated. 


nanoscale  patterns,12  and  influence  optical  properties  of  the 
inorganic  nanostructures. 2,10,11  For  example,  the  plasmon  reso¬ 
nance  frequency  of  metal  nanorods  depends  on  the  size  and 
shape  of  the  nanoparticles  as  well  as  on  the  presence  of 
molecular  recognition  groups,  such  as  peptides  or  surfactants 
attached  to  the  surface.  In  this  configuration,  binding  of  an 
external  agent  to  surface-attached  peptides  can  induce  sufficient 
shifts  in  UV/vis  absorption  to  support  sensor  capabilities.2 

A  major  challenge  lies  in  the  tremendous  number  of  potential 
peptides  and  lack  of  control  over  binding  affinities.  For  example, 
given  a  pool  of  20  natural  amino  acids  to  assemble  a  12-peptide, 
there  are  N  =  2012  =  1015  6  different  potential  peptides.  The 
screening  of  all  peptides  is  currently  impossible;  however,  phage 
display  techniques  allow  researchers  to  identify  peptide  se¬ 
quences  with  specific  binding  affinity  to  nanoparticles  and 
surfaces.  1,3,5~8,17_  9  this  procedure,  commercial  libraries  of 
several  billion  phages  with  distinct  peptide  sequences  expressed 
at  the  PHI  minor  coat  protein  were  screened  on  various  surfaces 
of  gold,  palladium,  and  silver  nanoparticles,  as  well  as  on 
surfaces  of  montmorillonite,  silica,  germanium  oxide,  tin  oxide, 
titania,  and  zinc  oxide.5,8  After  incubation  of  the  entire  library 
of  phages,  including  the  displayed  peptides  of  typically  8—12 
amino  acid  length,  with  the  nanoparticles  in  solution,  and  after 
several  washing  cycles  with  HEPES  buffer  solution  near  pH  = 
7,  the  most  strongly  binding  peptides  remain  bound  to  the 
inorganic  surfaces.  These  sequences  have  been  identified  by 
analysis  of  the  corresponding  phage  DNA  (some  shown  in  Table 
l).6  The  relative  binding  strength,  transmission  electron  mi¬ 
croscopy  (TEM),  IR  spectroscopy,  NMR  spectroscopy,  and 
thermochemical  data  were  recorded  to  characterize  the  peptide- 
modified  nanoparticles  in  comparison  to  the  peptides  in 
solution.1'13,17”19 

The  molecular  reasons  for  highly  specific  binding  in  the 
absence  of  thiol  groups  in  the  peptides,  however,  have  remained 
largely  a  matter  of  speculation.  Several  studies  attempted  an 
explanation  on  the  basis  of  experimental1,3'13,17'20  and  simula¬ 
tion  data;21'31  however,  evidence  about  the  acting  mechanisms 
is  often  vague,  and  modeling  studies  often  employed  excessively 
crude  approximations.  We  present  the  first  steps  using  molecular 
simulation  as  a  tool  to  understand  the  reasons  for  binding  of 
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1010),  Mineralization, 
Binding  Constants,  IR, 
NMR  experiments 


Figure  1.  Modeling  peptide— surface  interactions  and  biological  self- 
assembly  using  techniques  at  various  levels.  The  currently  distinguishable 
number  of  dodecapeptides  for  each  technique  is  indicated  in  parentheses. 
We  focus  on  classical  atomistic  simulation  in  comparison  with  experiment 
and  simulation  at  other  levels. 

selected  experimental  and  control  peptide  sequences,  including 
the  computation  of  adsorption  energies,  changes  in  peptide 
conformation  upon  adsorption,  and  the  strength  of  interaction 
of  amino  acid  residues  upon  interaction  with  the  surface  in 
aqueous  solution  (Figure  1).  In  contrast  to  earlier  work,  our 
approach  includes  an  accurate  representation  of  metal  surface 
energies32”36  using  a  Lennard-Jones  model  for  face-centered 
cubic  (fee)  metals  compatible  with  various  force  fields,29 
considers  the  dynamics  on  the  time  scale  of  nanoseconds  at  298 
K  as  opposed  to  static  calculations  at  0  K,  and  involves 
approximately  1000  water  molecules  as  opposed  to  unspecific 
implicit  water  models,  few  static  water  molecules,  or  no  water 
molecules  at  all.  Therefore,  temperature,  pH,  and  ionic  con¬ 
centration  are  tangible  quantities  in  the  simulation.  We  also 
benefit  from  an  efficient  new  method  for  the  computation  of 
binding  energies  for  the  peptides.  The  viability  of  this  approach 
is  supported  by  recent  improvements  in  the  accuracy  of  force 
field  parameters  for  inorganic  components29,37'41  and  the 
identification  of  polarization  effects  on  metal  surfaces42'46  as 
modest  contributions. 
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Therefore,  semiempirical  classical  simulation  in  close  relation 
to  experiment  and  other  theoretical  approaches  can  provide 
insight  into  the  nature  of  molecular  interactions  at  such 
inorganic— organic  interfaces  (Figure  1).  We  limit  the  investiga¬ 
tion  to  well-defined  even  {111}  and  {100}  surfaces;  the 
influence  of  the  shape  and  composition  of  the  nanoparticles  shall 
be  considered  elsewhere.  Further  background  information  on 
simulation  methods  can  be  found  in  section  SI  of  the  Supporting 
Information. 

The  paper  is  organized  as  follows.  In  section  2,  we  describe 
the  peptide  sequences  and  summarize  the  computational  ap¬ 
proach.  Sections  3  and  4  are  dedicated  to  the  presentation  of 
simulation  results  on  even  monometallic  surfaces  and  on 
bimetallic  surfaces.  Section  5  focuses  on  the  comparison  of 
results  to  experiment  and  explains  our  hypothesis  of  the 
specificity  and  strength  of  peptide  interactions  with  metal 
surfaces.  A  summary  is  presented  in  section  6.  The  Supporting 
Information  includes  a  description  of  the  computational  screening 
approach  (S2),  simulation  protocols  (S3),  and  additional  evidence 
on  the  binding  mechanism  from  earlier  simulations  (S4). 

2.  Peptide  Sequences  and  Computational  Details 

2.1.  Sequence  of  the  Peptides.  We  have  chosen  four  peptides 
named  A3,  Fig— Na3,  Pd2— Cl,  Pd4— Cl,  which  were  identified  as 
the  strongest  metal-binding  peptides  from  a  pool  of  ■-'lO9  peptides 
by  phage  display.7  8  The  neutral  peptide  A3  and  the  three-fold 
negatively  charged  peptide  Fig— Na3  at  pH  =  7  have  shown  affinity 
to  Au  and  Pd  nanoparticles.  The  peptides  Pd2— Cl  and  Pd4— Cl 
are  newly  reported  here  and  have  shown  affinity  to  Pd  nanoparticles. 
For  comparison,  we  also  include  the  peptide  Flgd— Na3,  deficient 
in  one  Asp  residue  compared  to  Fig— Na3,  as  well  as  the  arbitrarily 
chosen  peptides  GlylO  and  Pro  10  of  unknown  specificity  of 
adsorption  (Table  1).  We  assume  pH  =  7,  which  is  characterized 
by  protonation  or  deprotonation  of  the  corresponding  amino  acid 
side  chains  according  to  their  pA"a  values,  including  a  zwitterionic 
structure  of  the  peptides  with  a  protonated  N-terminal  end  and  a 
carboxylate  C-terminal  end.  For  amino  acids  whose  charges  do  not 
compensate  to  zero  within  the  peptide,  equivalent  amounts  of 
counterions  Cl-  or  Na+  are  added  to  maintain  charge  neutrality. 

2.2.  Efficient  Computation  of  Adsorption  Energies.  The 
quantitative  computation  of  adsorption  energies  requires  the  inclu¬ 
sion  of  explicit  water  molecules  and  ionic  solutes  due  to  the 
formation  of  spatially  oriented  hydrogen  bonds  which  change  their 
orientation  when  peptides  (or  other  adsorbates)  move  from  the 
solution  to  the  surface.  We  utilize  a  new,  efficient  approach  which 
involves  the  simulation  of  four  systems,  (1)  the  surface-solute— 
solvent  system  of  average  energy  E\,  (2)  the  solute— solvent  system 
of  average  energy  £2,  (3)  the  solvent  system  of  average  energy  £3, 
and  (4)  the  surface-solvent  system  of  energy  £4  under  consideration 
of  equal  molecular  volumes  for  each  component  under  NVT 
conditions.  The  adsorption  energy  is  calculated  from  the  average 
energy  of  each  of  the  four  systems  as  Eads  =  E\  -  E2  +  £3  -  £4. 
Each  simulation  requires  an  accurate  summation  of  Coulomb 
interactions  on  the  order  of  <0.1  kcal/mol  uncertainty  for  the  entire 
system  as  well  as  precise  temperature  control  or  energy  corrections 
to  a  reference  temperature,  respectively.  Details  can  be  found  in 
section  S2  of  the  Supporting  Information. 

2.3.  Models,  Force  Field  Parameters,  and  Simulation  Pro¬ 
tocol.  Fully  atomistic  models  of  the  peptides,  metal  surfaces,  and 
water  were  prepared  using  the  Hyperchem47  and  Materials  Studio 
Visualizers.48  Box  dimensions  were  typically  3  nm  x  3  nm  x  z 
whereby  the  additive  box  height  z  amounts  to  >1.2  nm  for  metal 

(46)  Guidelli,  R.;  Schmickler,  W.  Electrochim.  Acta  2000,  45,  2317-2338. 

(47)  HyperChem  7.5',  Hypercube  Inc.:  Gainesville,  FL,  2006. 

(48)  Materials  Studio  4.0,  Discover  Program,  Cerius2,  and  User  Guide; 

Accelrys,  Inc.:  San  Diego,  CA,  2005. 


A3 


Flg-Na3 


Flg-Na2 


Figure  2.  Representative  snapshots  of  the  peptides  in  solution.  The  peptides 
A3,  Fig— Na3,  and  Flgd— Na2  assume  unspecific  3D  conformations,  and  a 
tendency  toward  intramolecular  interactions  between  the  charged  —  CCF- 
groups  and  H3N+— /Na+  ions  is  seen  for  Fig— Na3  and  Flgd— Na3.  The 
peptide  GlylO  prefers  an  extended  conformation  and  is  conformationally 
very  flexible.  Polyproline  ProlO  assumes  a  helical  structure. 

substrates,  and  to  >3.5  nm  for  water  and  aqueous  solutions  of  the 
peptides.  The  size  was  chosen  such  that  25%  to  40%  of  the  surface 
area  is  covered  by  the  peptide,  close  to  experimental  conditions. 
The  consistent  valence  force  field  (CVFF)48'49  extended  for  accurate 
Lennard-Jones  parameters  for  fee  metals29  was  employed  in  the 
simulation  of  the  metal-water— peptide  hybrid  systems.  The  pa¬ 
rameters  for  fee  metals  reproduce  surface  energies  and  interface 
energies  in  quantitative  agreement  with  experiment  (1  —  10%),  down 
from  deviations  up  to  500%  in  earlier  force  fields.  The  Discover 
program48  as  well  as  LAMMPS50  were  employed  to  carry  out 
molecular  dynamics  simulations.  Models,  force  field,  simulation 
protocols,  and  accuracy  are  described  in  section  S3  of  the 
Supporting  Information. 

3.  Adsorption  of  Peptides  on  Even  {111}  and  {100} 
Surfaces  of  Au  and  Pd 

The  peptides  were  relaxed  in  aqueous  solution  (Figure  2)  and 
on  the  metal  {111}  and  { 100}  surfaces  (Figure  3)  by  molecular 
dynamics  simulation.  Corresponding  adsorption  energies  (Table 
2)  indicate  that  peptides  are  generally  attracted  to  the  metal 
surfaces,  supported  by  the  observation  of  direct  contact  with 
the  metal  surface  or  the  presence  of  one  or  two  water  interlayers 
in  the  models  (Figure  3).  The  order  of  magnitude  of  computed 
binding  energies  for  the  oligopeptides,  between  0  and  —100 
kcal/mol,  is  in  the  range  of  noncovalent  interactions  (<100  kcal/ 
mol)  and  often  stronger  than  that  of  a  simple  collection  of 
hydrogen  bonds  (1—4  kcal/mol).  The  adsorption  strength 
correlates  with  the  amount  of  direct  contact  with  the  metal 
surface  and  differs  significantly  between  the  {111}  and  the 
{ 100}  surfaces.  In  the  peptides,  we  can  identify  binding  residues 
versus  nonbinding  residues.  Strongly  binding  peptides  contain 
several  binding  residues,  while  the  remaining  residues  are  less- 
binding  or  nonbinding.  In  Flgd— Na2,  the  elimination  of  a  single 
amino  acid  (Asp)  in  the  peptide  sequence  compared  to 
Flgd— Na3  leads  to  a  modest  reduction  in  binding  strength  (Table 
2). 


(49)  Dauber-Osguthorpe,  P.;  Roberts,  V.  A.;  Osguthorpe,  D.  J.;  Wolff,  J.: 
Genest,  M.;  Hagler,  A.  T.  Proteins:  Struct.,  Fund.  Genetics  1988,  4, 
31-47. 

(50)  Plimpton,  S.  J.  J.  Comp.  Phys.  1995,  117,  1-19. 
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Figure  3.  Representative  snapshots  of  the  peptides  A3,  Fig— Na3,  Flgd— Na2,  GlylO,  and  ProlO  on  the  metal  surfaces,  viewed  along  different  perspectives 
to  support  a  spatial  impression.  Some  counterions  are  not  shown  due  to  diffusion  out  of  the  displayed  region. 


3.1.  Direct  Contact  to  the  Metal  Surface  and  Adsorption 
Strength.  We  define  direct  contact  as  a  distance  of  the  atoms 
of  residues  or  chain  fragments  <3  A  from  the  metal  surface 
atomic  plane  for  the  majority  of  the  time,  with  typical  distances 
of  2.5  A.  In  cases  of  strong  binding,  individual  side  groups  or 
backbone  residues  are  found  in  direct  contact  with  the  metal 
surface,  without  intermediate  water  molecules  between  the 
peptide  and  the  metal  surface  (Figure  3).  In  cases  of  weak 
binding,  there  is  no  such  direct  contact,  and  multiple  “layers” 
of  water  molecules  accumulate  between  the  metal  surface  and 
the  peptide.  For  a  single  water  layer  in  between,  the  atoms  of 
the  residues  or  chain  fragments  nearest  to  the  metal  surface 
atomic  plane  are  typically  found  within  a  distance  of  5  A  for 
the  majority  of  the  time.  For  multiple  water  layers  in  between, 
the  average  distance  is  accordingly  greater. 

For  all  peptides,  adsorption  is  stronger  on  the  {111}  surface 
compared  to  the  { 100}  surface.  The  adsorption  energy  on  { 100} 
surfaces  amounts  only  to  ~10%  compared  to  {111}  surfaces 
for  the  peptides  A3  and  Fig— Na3.  The  control  peptide  GlylO 


also  shows  a  strong  adsorption  on  {111}  surfaces  and  high 
selectivity.  The  control  peptide  ProlO  is  the  least  binding  peptide 
in  all  instances,  although  moderate  attraction  is  also  seen  to 
{111}  surfaces.  The  visualization  supports  that  oligoprolines 
may  be  suitable  to  form  loops  to  connect  more  specialized 
peptide  fragments  (Figures  2  and  3).  In  binding  peptides,  at  least 
a  few  binding  amino  acids,  groups  of  amino  acids,  or  the  entire 
peptide  are  in  direct  contact  with  the  metal  surface  for  most  of 
the  time,  which  accounts  for  the  majority  of  the  binding  energy. 
The  correlation  between  independent  conformations  of  the  same 
peptide  and  corresponding  adsorption  energies  further  indicates 
cooperativity  between  adjacent  amino  acids:  the  presence  of  less- 
binding  residues  near  binding  residues  in  the  peptide  chain  as 
well  as  the  polymerization  of  individually  less  strongly  binding 
residues  to  homopeptides  (GlylO,  ProlO)  can  increase  the  extent 
of  direct  contact  with  the  surface  and  the  overall  binding  energy. 

3.2.  Amino-Acid-Specific  Contributions  to  Adsorption 
Strength.  Analysis  of  the  MD  trajectories  provides  the  first 
amino-acid-specific  trends  to  guide  in  the  design  of  peptides 
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Table  2.  Adsorption  Energies  of  the  Peptides,  in  kcal/mol  (±5  kcal/mol  or  up  to  ±10%), 
with  the  Surface3 

and  a  Compilation  of  Residues  in  Direct  Contact 

peptide 

metal  surface 

A3 

Fig — Na3 

Flgd— Na2 

GlylO 

ProlO 

Au  {100} 

-9 

2Tyr,  4Ser,  8Pro  (partly), 
9Met  (main  chain), 
loPro 

0 

2  C02-  of  Asp 

—6 

2  C(V  of  Asp  and  7Lys 

-4 

none 

0 

none 

Au  { 1 1 1 } 

-63 

first  five  residues 
iAla— 2Tyr— 3Ser— 4 

Ser— 5Gly,  gPro 
(partial),  nPro,  and 
nPhe  (incl.  CC>2_) 

-63 

all  residues;  2Tyr,  1—2 
C02_  of  i)4_7Asp,  8Lys 
partly  detached;  Na+ 

— CO2-  bridges 

-45 

iAsp  and  3Lys  to  7Lys; 
NH3+  of  iAsp,  3>7Lys 
and  1—2  CC>2_  of  Asp 
partly  detached 

-60 

6—7  consecutive  Gly 
residues 

-23 

2—3  Pro  and  additional 
partial  contacts 

Pd  {100} 

-2 

iAla,  2Tyr,  4Ser,  8Pro 
(partly),  gMet  (main 
chain),  and  ioPro 

+3 

1—2  CC>2~  of  Asp  and 
3Lys 

+4 

NH3+  of  iAsp,  C02-  on 
7Lys,  sometimes  CC>2_ 
of  Asp 

-27 

1—2  Gly  residues 

-13 

1  Pro  (partial)  and  C02- 
of  ioPro  (partial) 

Pd  { 1 1 1 } 

-77 

first  five  residues 
iAla— 2Tyr— 3Ser— 4 

Ser— 5Gly,  gPro 
(partial),  nPro,  and 
i2Phe 

-85 

all  residues;  NH3+  of 
iAsp,  3,8Lys,  and  1—2 
C02-  (Asp)  less 
attached,  Na+— CC>2_ 
bridges 

-60 

all  residues  except  3Lys; 

1  C02-  of  Asp  and 
NH3+  of  iAsp  partly 
detached 

-75 

7—9  residues  incl.  C02 
of  ioGly 

-40 

5—6  Pro  and  additional 
partial  contacts 

The  subscript  to  the  left  indicates  amino  acid  positions  within  the  peptide. 


Figure  4.  Top  view  of  representative  snapshots  of  (a)  Fig— Na3  and  (b)  Flgd— Na2  on  a  Pd  {111}  surface,  (a)  In  Fig— Na3,  all  residues  are  in  direct  contact 
with  the  surface,  and  ions  pairs  (1  and  2)  are  highlighted,  (b)  Flgd— Na2  contains  one  Asp  residue  less  and  forms  different  internal  ion  pairs  (1—3).  The  side 
chain  of  3Lys  in  feature  3  is  detached  from  the  surface. 


(Table  2). 51  In  A3,  the  phenyl  groups  in  2Tyr  and  ,2Phe  attach 
preferably  to  the  planar  {111}  Pd  surface,  interconnected  by 
less  tightly  bound  Ala,  Pro,  and  Met  residues.  On  the  Au  {111} 
surface,  the  phenyl  ring  of  i2Phe  is  less  tightly  bound  compared 
to  that  on  the  Pd  {111}  surface.  The  exact  analysis  of  the 
adsorption  energy  of  individual  residues  is  challenging  due  to 
changes  in  solvation  and  conformation  upon  contact  with  the 
surface.  Nevertheless,  the  two  aromatic  units  and  their  sur¬ 
rounding  region  in  the  peptide  chain  can  be  regarded  as 
anchoring  residues  and  account  for  much  of  the  binding  energy. 
This  view  concurs  with  a  less  intense  adsorption  of  A3  on  the 
Au  {111}  surface  compared  to  the  Pd  {111}  surface.  Direct 
contact  with  fewer  residues  of  A3  is  observed  on  the  Pd  {100} 
and  Au  {100}  surfaces  (Table  2).  The  adsorption  strength  of 
the  individual  amino  acids  is  also  significantly  weaker. 

In  Fig— Na3,  all  residues  of  the  peptide  are  in  direct  contact 
with  the  Pd  {111}  surface,  mostly  driven  by  the  2Tyr  residue 
as  well  as  3—4  C02~  groups  of  14_7Asp  and  the  C-terminal 
sLys  (Figure  4a).  The  3,8Lys  side  chains  and  1—2  C02  groups 


(51)  A  systematic  study  of  the  adsorption  of  single  amino  acids  on  a  metal 
surface  in  solution  will  be  reported  separately. 


are  not  in  direct  contact  for  part  of  the  time.  Notable  are 
intramolecular  bridges  between  Na+  ions  and  C02  groups 
(features  1  and  2).  On  the  Auflll}  surface,  the  adsorption 
geometry  is  similar  except  for  somewhat  reduced  contact  with 
2Tyr  and  lower  adsorption  energy.  On  the  Pd  and  Au  {100} 
surfaces,  one  or  two  C02  groups  of  the  Asp  residues  or  3Lys 
are  in  contact  with  the  metal  surface,  yet  essentially  no 
adsorption  results. 

In  Flgd— Na2,  all  groups  except  3Lys  directly  interact  with 
the  Pd  {111}  surface  (Figure  4b).  Complexes  between  a  sodium 
ion  with  C02  groups  can  be  seen  (feature  1);  the  other  sodium 
ion  floats  at  a  distance  in  solution.  7Lys  coordinates  with  5Asp 
(feature  2),  and  3Lys  coordinates  with  4Asp  (feature  3).  On  the 
Au  {111}  surface,  2Tyr  is  not  in  contact  with  the  surface,  while 
the  rest  of  the  peptide  including  3Lys  is  bound.  The  binding 
strength  is  mainly  associated  with  C02  groups  in  Asp  and  7Lys 
and  sometimes  another  C02  group  of  6Asp.  The  adsorption  is 
weaker  compared  to  that  on  Pd  {111},  even  though  the  amount 
of  direct  contact  is  similar.  Less  contact  is  seen  on  the  Au  and 
Pd  {100}  surfaces  (Table  2),  and  the  adsorption  is  small  to 
negligible. 
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The  reference  peptides  GlylO  and  Pro  10  also  exhibit  interest¬ 
ing  trends.  Most  notably,  the  binding  potential  of  GlylO  is  high. 
Seven  to  nine  Gly  residues,  including  the  C-terminal  CO? 
group,  bind  to  Pd  {111},  and  some  residues  including  6Gly  near 
the  center  of  the  chain  favor  a  water  interlayer  and  detach  from 
the  surface  (Figure  3).  Six  or  seven  consecutive  Gly  residues 
toward  the  C-terminal  end  bind  to  Au  {111},  and  the  adsorption 
energy  is  close  to  that  of  A3  and  Fig— Na3.  Also,  the  Pd  { 100} 
surface  attracts  a  few  Gly  residues,  while  Au  {100}  shows  no 
adsorption. 

Pro  10  is  essentially  inert  against  the  Au  {100}  surface  and 
shows  some  adsorption  through  two  or  three  Pro  residues  on 
the  Au  {111}  surface.  It  interacts  slightly  with  Pd  {100},  which 
results  in  minor  adsorption,  and  substantially  with  Pd  {111} 
through  five  or  six  residues  and  other  partial  contacts,  including 
the  C-terminal  CO?  group.  Flowever,  the  proline  helix  resists 
flattening  so  that  the  adsorption  energy  remains  less  than  half 
the  value  compared  to  other  peptides. 

3.3.  Differences  in  Secondary  Structure  upon  Adsorption. 
Adsorption  of  the  peptides  is  also  related  to  differences  in 
solvation  in  the  presence  of  the  surface  compared  to  the  state 
in  solution.  The  peptides  are  short  but  of  a  critical  length  to 
observe  secondary  structure  formation. 

Peptide  A3  in  solution  is  of  a  nonspecific,  flexible  conforma¬ 
tion  which  essentially  converts  into  a  flat-on  (2D)  structure  in 
the  presence  of  {111}  surfaces.  A  few  residues  close  to  Pro  in 
the  middle  section  of  the  peptide  are  not  directly  bound  and 
retain  flexibility.  On  the  {100}  surfaces,  the  peptide  exhibits  a 
3D  structure  closer  to  the  solution  state.  Peptide  Fig— Na3  in 
solution  forms  intramolecular  —  C02”"'H3N+—  bridges  be¬ 
tween  sLys  and  s,4Asp  as  well  as  between  3Lys  and  6jAsp,  which 
lead  to  a  3D  shaped  (and  concentration-dependent)  conformation 
(Figure  2).  On  the  {111}  surfaces,  these  intramolecular  bridges 
are  given  up  in  favor  of  a  flat-on,  2D  shape  which  involves 
several  C02  •••Nal  bridges  (Figure  4a).  On  the  less  attractive 
{100}  surfaces,  more  3D  features  are  preserved.  Peptide 
Flgd— Na2,  missing  one  Asp  residue  compared  to  Fig— Na3, 
exhibits  less  intramolecular  bridging  in  solution.  Upon  adsorp¬ 
tion  on  the  {111}  surfaces,  it  approaches  a  2D  backbone  except 
for  one  or  two  residues,  supported  by  stabilizing  ion  pairs 
(Figure  4b).  The  solution  structure  is  also  preserved  to  a  higher 
extent  on  the  {100}  surfaces.  Peptide  GlylO  in  solution  prefers 
an  extended  conformation  with  significant  conformational 
flexibility  (Figure  2).  Upon  adsorption,  the  peptide  appears 
mostly  “glued”  to  the  Pd  {111}  surface,  less  so  to  the  Au  {111} 
surface,  and  maintains  significant  conformational  flexibility  on 
the  {100}  surfaces  (Figure  3).  Peptide  ProlO  in  solution  forms 
a  mixture  of  polyproline  II  and  polyproline  I  helices  due  to  the 
steric  demand  of  the  pyrrolidine  rings  (Figure  2).  Upon 
adsorption,  the  conformation  is  biased  toward  contact  with  the 
surface.  Particularly  on  Pd  {111},  more  than  half  the  residues 
are  in  contact  with  the  surface,  less  so  on  Au  {111}.  Minor 
contacts  are  seen  on  the  Pd  {100}  surface  and  no  contact  with 
the  Au  {100}  surface  (Table  2). 

Significant  conformational  changes  such  as  a  bias  toward  2D 
structures  and  greatly  reduced  mobility  of  individual  groups  on 
{111}  surfaces  upon  adsorption  indicate  a  significant  loss  in 
entropy.  While  it  is  time-consuming  to  quantify  the  entropy 
balance  computationally,  a  crude  estimate  can  be  based  on  loss 
in  rotational  degrees  of  freedom  of  the  peptide  and  an  increase 
in  lateral  mobility  of  bound  water  such  as  during  a  melting 
transition.  The  assumption  of  four  rotatable  bonds  with  three 
equal  rotational  states  each  per  amino  acid  (main  chain  + 


Table  3.  Data  on  the  Polarity  of  Pd-Au  Interfaces  for  Use  in  an 
Extended  Born  Model,38  in  MJ/mol 


property 

Pd 

Au 

sublimation  energy" 

0.38 

0.37 

first  ionization  energy" 

0.80 

0.89 

second  ionization  energy" 

1.88 

1.98 

first  electron  affinity " 

0.054 

0.223 

enthalpy  of  alloy  formation*7 

~ 

0.02  <AH<0 

“  Reference  35.  h  References  52—54. 


average  side  chain)  for  a  10-peptide  results  in  a  conformational 
entropy  of  S  =  10k  In  34,  equal  to  a  maximum  loss  in  entropic 
freedom  of  —  TAS  ~  40 kT  «  24  kcal/mol.  The  concomitant 
release  of  25  water  molecules,  which  approximately  corresponds 
to  the  surface  area  covered  by  an  extended  10-peptide,  results 
in  a  maximum  entropy  increase  of  25  AS(,  where  ASf  corresponds 
to  the  fusion  entropy35  and  leads  to  —SAT  =  —  39  kcal/mol. 
We  observe  in  the  simulation  that  certainly  not  all  degrees  of 
freedom  are  lost,  particularly  so  for  bound  water  molecules,  so 
that  actual  values  may  be  less  than  half  of  the  maximum  and 
significantly  compensate  each  other.  Therefore,  computed 
binding  energies  in  Table  2  are  likely  acceptable  approximations 
of  binding  free  energies  as  well. 

4.  Adsorption  of  Peptides  on  a  Bimetallic 
Pd{1 1 1 } — Au{1 11}  Interface 

We  further  analyzed  the  adsorption  of  the  peptides  A3, 
Fig— Na3,  Flgd— Na2,  Pd2Cl,  and  Pd4Cl  on  even  bimetallic 
Pd— Au  {111}  surfaces  with  a  smooth,  idealized  boundary 
between  Pd  and  Au.  The  boundary  creates  an  electrostatic 
potential  difference  between  the  two  metals,  and  the  method  of 
Heinz  and  Suter38  was  employed  for  an  estimate  of  the  polarity, 
since  quantitative  experimental  data  were  not  available. 

4.1.  Polarity  of  the  Bimetal  Interface.  The  ground-state 
electron  configurations  of  Pd  and  Au  are  [Pd]  5s°4d10  and  [Au] 
bsHf'^d10,  in  which  the  partially  occupied  6s  orbital  of  Au 
indicates  its  action  as  an  electron  acceptor.  Atomization  and 
ionization  energies  confirm  this  trend  and  permit  an  estimate 
of  the  atomic  charges  at  the  interface  using  an  extended  Born 
model  (Table  3). 38  The  enthalpy  of  alloy  formation,  and 
therefore  the  enthalpy  of  formation  of  the  Pd— Au  interface,  is 
close  to  zero.52”54  The  analysis  for  different  atomic  charges 
(±0e,  ±0.5e,  ±1.0e,  ±2.0e)  in  a  bulk  mixture  and  for  a  single 
flat  interface  suggests  a  polarity  near  ±le  in  a  3D  bulk  alloy 
and  near  ±0.3e  at  a  2D  flat  interface  (Figure  5). 

While  we  note  a  significant  difference  in  atomic  charges 
between  a  bulk  mixture  and  a  smooth  interface,  we  employ 
charges  of  +0.3e  for  Pd  and  — 0.3e  for  Au  atoms  (±0.1e)  at 
the  flat  interface  (all  other  metal  atoms  ±0.0e).  The  assignment 
agrees  with  X-ray  measurements  of  the  electron  deformation 
density  of  reference  structures  of  the  form  A^.38 

4.2.  Trends  in  Adsorption  Energies.  We  computed  the 
adsorption  energy  of  the  peptides  on  the  Pd— Au  {111} 
bimetallic  surface  (Table  4).  The  presence  of  distinct  surface 
areas  of  Au  and  Pd  allowed  monitoring  the  preferred  location 
of  the  peptides  in  the  course  of  the  simulation  (Figure  6). 
Different  start  positions  were  chosen  arbitrarily  between  the  Au 
and  Pd  surface  which  do  not  affect  the  final  result  after  sufficient 


(52)  Tayler,  J.  B.  Proc.  Phys.  Soc.  London  1899,  17,  194-201. 

(53)  Niessen,  A.  K.;  de  Boer,  F.  R.;  Boom,  R.;  de  Chatel,  P.  F.;  Mattens, 
W.  C.  M.  CALPHAD  1983,  7,  51-70. 

(54)  Meschel,  S.  V.;  Kleppa,  O.  J.  J.  Alloys  Compd.  2003,  350,  205-212. 
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Figure  5.  Estimate  of  atomic  charges  at  Pd— Au  interfaces.  Compared  to  (a)  a  3D  homogeneous  distribution,  (b)  a  flat  interface  leads  to  less  efficient  charge 
transfer.  The  gain  in  Coulomb  energy  per  charged  atom  amounts  only  to  1/8  for  the  same  atomic  charges. 


Table  4.  Adsorption  Energies  at  the  Pd-Au  {111}  Bimetal  Surface,  in  kcal/mol  per  Peptide  (=1=5  kcal/mol  or  up  to  ±10%),  and  a  Compilation 
of  Residues  in  Direct  Contact  with  the  Surface 


metal  surface 

peptide 

Pd-Au  (111) 

Pd  (111) 

Au  {111} 

A3 

-86 

most  of  jAla  to  6Ala,  gPro  and  ioPro  (partly),  nPro  and  i2Phe 

-77 

-63 

Fig— Na3 

-74 

all  residues,  except  sometimes  NH3+  of  iAsp,  side  chains  in  1—2  Asp 
and  in  3;sLys 

-85 

-63 

Flgd— Na2 

Pd2— Cl 

Pd4— Cl 

-63 

iAsp,  2Tyr,  4Asp  to  7Lys,  side  chains  in  1—2  residues  (Asp,  7Lys) 

sometimes  detached,  ion  pair  between  3Lys  and  6Asp 

-78 

all  residues;  iAsn  inch  NH3+,  3Met,  5Leu,  9Gly,  and  ^His  inch  CO2- 

sometimes  detached 

-93 

all  residues  except  7Pro;  sVal,  gThr,  gLeu,  and  i2Leu  sometimes  detached 

-60 

-45 

simulation  time.  In  some  cases,  lateral  offset  of  the  peptide  and 
continued  simulation  on  the  order  of  several  nanoseconds  were 
required  to  identify  energetically  preferred  locations.  All  peptides 
remain  spatially  flexible  and  gravitate  toward  the  Pd  domain  of 
the  surface.  Due  to  polarity,  the  Pd— Au  bimetal  interface  acts 
as  a  preferred  adsorption  site  for  A3  and  lowers  the  (negative) 
adsorption  energy  further  by  about  10  kcal/mol  in  comparison 
to  the  most  attractive  monometallic  {111}  surface.  However, 
less  attraction  is  observed  for  charged  Fig— Na3  (Table  4).  In 
every  case,  adjustments  in  conformation  of  the  peptide  backbone 
occur  to  adapt  to  the  polarity  of  the  Pd— Au  interfacial  region. 
Patterns  of  binding  and  less-binding  residues  are  similar  to  those 
observed  on  monometallic  surfaces,  and  the  adsorption  strength 
is  of  similar  magnitude.  In  Pd2— Cl  and  Pd4— Cl,  Phe,  His,  Arg, 
and  Asn  are  the  most  strongly  binding  residues. 

Peptide  A3  translates  to  the  Pd— Au  interface  in  the  initial 
stages  of  the  simulation  and  remains  there  for  the  rest  of  the 
time.  Almost  the  entire  peptide  is  in  direct  contact  with  the  metal 
surface,  and  the  adsorption  strength  is  increased  relative  to  that 
observed  on  Pd  {111}  and  Au  {111}  surfaces.  The  orientation 
of  polar  functional  groups  in  the  backbone  and  in  the  side  chains 
complementary  to  the  polar  Pd+0-3  •  •  •  Au~03  interface  may 
account  for  this  effect,  for  example,  — NH— ,  —OH  groups  of 
2Tyr  and  3Ser  near  interfacial  Au,  and  carbonyl  oxygen  atoms 
of  io,i iPr°  near  interfacial  Pd  (Figure  6).  Peptide  Fig— Na3  also 
prefers  the  Pd  domain,  and  all  residues  assume  direct  contact 
with  the  metal  similar  to  the  observation  on  the  Pd  {111} 
surface.  However,  a  reduction  in  adsorption  strength  by  ca.  10 
kcal/mol  was  found.  For  example,  C02  groups  of  i,6Asp  and 
8Lys  are  found  near  interfacial  Pd  (+0.3e),  and  — NH3+  and 
—OH  of  3Lys  and  2Tyr  are  found  in  the  vicinity  of  interfacial 
Au  (— 0.3e).  The  peptide  Flgd— Na2  is  also  attracted  to  the  Pd 
domain  of  the  surface,  and  similar  characteristic  orientation  of 
polar  groups  along  the  Pd— Au  interface  is  seen  (Figure  6). 

The  peptide  Pd2— Cl  exhibits  an  extended  conformation  and 
no  intramolecular  ion  pairs  in  solution.  On  the  Pd— Au  bimetal 
{111}  surface,  the  entire  peptide  is  in  direct  contact  with  the 
surface  and  adsorbs  strongly  with  —78  kcal/mol.  Simulation 
with  different  start  structures  indicates  that  the  residues  iAsn, 


2Phe,  7Arg,  and  io,i2His  contribute  particularly  to  adsorption.  A 
preference  for  the  Pd  domain  of  the  bimetal  surface  is  noticed, 
and  the  polarity  of  the  Pd— Au  interface  supports  the  adsorption 
of  some  residues  in  its  vicinity  (Figure  6).  The  peptide  Pd4— Cl 
is  also  of  an  extended  conformation  in  solution  without 
intramolecular  ion  pairs.  Upon  adsorption,  direct  contact  with 
the  bimetal  is  made  through  all  residues  but  7Pro.  Hydrophobic 
residues  are  partly  detached  from  the  surface  (Table  4),  and 
the  adsorption  energy  is  very  high  at  —93  kcal/mol,  even  though 
the  peptide  extends  across  Pd  and  Au  domains.  Simulations  with 
different  start  structures  indicate  that  3Asn,  6,nHis,  and  ioArg, 
as  well  as  i,8Thr,  contribute  particularly  to  adsorption. 

5.  Comparison  with  Experiment  and  Binding 
Mechanism 

The  simulation  provides  semiquantitative  measures  of  adsorp¬ 
tion  energies  and  insight  into  the  rich  interplay  of  factors  that 
contribute  to  binding.  Nevertheless,  we  need  to  consider 
limitations  of  the  model.  We  neglect  contributions  due  to 
polarization  and  charge  transfer  at  the  metal  interface  (covalent 
bonding  contributions),  finite-size  effects  owed  to  limitations 
in  box  size  to  3—5  nm  and  application  of  periodic  boundary 
conditions,  shortcomings  of  the  force  field  models,  and  the 
existence  of  other  peptide  conformations  which  may  not  have 
been  captured  due  to  long  relaxation  times  on  the  metal  surfaces. 
Yet  the  body  of  results  is  consistent  with  numerous  experimental 
data  and  permits  substantial  conclusions  on  the  mechanism  of 
binding  as  discussed  in  the  following. 

5.1.  Experimental  Evidence  for  Adsorption  of  Peptides  on 
Pd,  Au,  and  Pd— Au  Bimetal  Surfaces.  Brown  et  al.1  identified 
noncovalent  adhesion  of  peptides  to  Au  {111}  faces  by  protein- 
induced  growth  of  Au  nanocrystals.  Fu  et  al.4  observed  the 
adsorption  of  Pd  and  Au  nanoparticles  of  4—5  nm  average  size 
onto  a  fibrillar  dodecapeptide  (Tl,  Arg-Gly-Tyr-Phe-Trp-Ala- 
Gly-Asp-Tyr-Asn-Tyr-Phe).  The  metal  particles  were  obtained 
by  reduction  of  PdCl2  and  HAuCI4  in  solution,  which  typically 
leads  to  {111}  surfaces.  The  high  content  of  aromatic  residues 
Tyr,  Phe,  and  Trp  in  peptide  Tl  was  associated  with  strong 
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Figure  6.  Representative  snapshots  of  peptides  on  a  Pd— Au  {111)  bimetal 
surface.  A  preference  for  Pd  domains,  the  orientation  of  polar  functional 
groups  along  the  Pd+0  3  •  •  •  Au  0  3  interface,  and  the  position  of  residues  on 
the  hexagonal  surface  can  be  seen.  Water  molecules  are  partially  shown  to 
indicate  the  presence  of  solvent  and  maintain  a  clear  view.  The  motion  of 
larger  peptides  is  restrained  by  finite  box  size. 

adsorption  of  metal  nanoparticles,  and  metallic  double-helical 
structures  were  obtained.  This  observation  coincides  well  with 
the  strong  attraction  of  aromatic  side  chains  to  both  Au  {111} 
and  Pd  {111}  surfaces  in  the  simulation  (Tables  2  and  4).  Fu  et 
al.  also  mention  the  common  adsorption  of  anionic  groups  such 
as  CCL  on  colloidal  metal  nanoparticles.  This  observation 


concurs  with  the  binding  of  Asp  side  chains  in  Fig— Na3  and 
Flgd— Nai  to  Au  {111}  and  Pd  {111}  surfaces,  and  possibly 
the  proximity  of  the  Asp  side  chains  to  Au  {100}  and  Pd  {100} 
surfaces  in  the  simulation  (Tables  2  and  4). 

However,  measurements  of  the  adsorption  energy  of  indi¬ 
vidual  amino  acids  and  peptides  have  yet  rarely  been  reported. 
Li  et  al.  measured  the  adsorption  energy  of  pure  acetic  acid  on 
Au— Pd  {111}  alloy  surfaces  by  temperature-programmed 
desorption  and  reflection— absorption  infrared  spectroscopy.20 
While  a  quantitative  comparison  to  aqueous  solution  is  not 
possible,  acetic  is  similar  to  Asp  side  chains  in  the  Fig— Na3 
and  Flgd— Na3  peptides,  or  to  Gly  in  protic  environment.  When 
the  mole  fraction  of  Au  is  greater  than  0.5  in  the  bimetal,  acetic 
acid  adsorbs  and  desorbs  reversibly  with  ~12  kcal/mol.  When 
the  Au  fraction  is  below  0.33,  acetic  acid  coordinates  as  an  rj2 
ligand  at  room  temperature  or  decomposes  on  the  surface.20  The 
measured  adsorption  energy  is  qualitatively  consistent  with 
computed  adsorption  energies  for  Fig— Na3  and  Flgd— Na2 
peptides  in  water  solution  on  the  Pd— Au  {111}  surface  (Table 
4):  approximately  4—5  CC>2~  groups  coordinate  to  the  surface, 
which  would  amount  to  48—60  kcal/mol  according  to  Li  et  al.; 
additional  adsorption  strength  is  gained  through  the  Tyr  residue 
and  the  peptide  backbone  which  accounts  for  the  total  computed 
adsorption  energy  of  45—85  kcal/mol  (Table  4).  We  have  also 
seen  a  tendency  toward  rf2  coordination  of  the  Asp  side  chains 
on  Pd  {111}  surfaces  in  the  simulation  (Figures  4  and  6)  in 
agreement  with  observations  for  pure  acetic  acid.20 

The  interaction  of  peptides  with  Ag,  Au,7,8  and  Pd8  salt 
solutions  was  investigated  by  Slocik  et  al.  A3  and  Fig— Na3 
reduce  HAuC’L,  possibly  through  oxidation  of  the  2Tyr  residue, 
and  attach  to  the  surface  of  the  emerging  metal  nanoparticle 
with  thermodynamically  stable  {111}  faces32-34'55  under  non- 
oxidative  conditions.56'57  On  the  surface  of  the  Au  nanoparticles, 
A3  leads  to  a  nanoparticle  solution  and  Fig— Na3  to  nanoparticle 
precipitation.  The  average  size  (~diameter)  of  the  fee  nanoc¬ 
rystalline  Au  particles  of  ~13  nm  resulted  in  a  narrow  size 
distribution  when  the  surface  was  covered  with  A3  and  a  wider 
size  distribution  of  irregularly  shaped  particles  when  covered 
with  Fig— Na3.  Simulation  shows  a  clear  preference  for  adsorp¬ 
tion  onto  Au  {111}  surfaces  over  Au  {100}  surfaces,  with 
adsorption  energies  of  —63  and  —63  kcal/mol  compared  to  —9 
and  0  kcal/mol  for  A3  and  Fig— Na3,  respectively  (Table  2). 
The  peptides  of  2.0— 3.5  nm  length  are  also  small  in  relation  to 
~13  nm  diameter  of  the  grown  particles,  so  an  assumption  of 
even  surfaces  (neglect  of  curvature)  in  the  simulation  may  be 
acceptable.  Separately,  the  peptides  were  incubated  with  Au 
nanoparticles  of  30  nm  diameter,  and  fluorescence  measurements 
of  2Tyr  (emission  303  nm)  indicated  a  stronger  relative 
adsorption  of  A3  vs  Fig— Na3  of  6.8  vs  5.1  with  an  uncertainty 
of  0.4  unit.7  The  simulation  indicates  equal  values;  however, 
adsorption  of  the  longer  A3  peptide  could  be  slightly  underes¬ 
timated  due  finite  size  effects  in  the  comparatively  small 
simulation  box.  The  difference  between  dispersion  of  A3-derived 
nanoparticles  and  precipitation  of  Fig— Na3-derived  nanoparticles 
appears  to  be  related  to  different  charge  states,  in  addition  to 
differences  in  adsorption  energy:  A3-coated  nanoparticles  are 
neutral  and  Fig— Na3  coated  nanoparticles  are  significantly 


(55)  Buffat,  P.-A.;  Flueli,  M.;  Spycher,  R.;  Stadelmann,  P.;  Borel,  J.-P. 
Faraday  Discuss.  1991,  92,  173-187. 

(56)  Wiley,  B.;  Sun,  Y.-G.;  Chen,  J.-Y.;  Cang,  H.;  Li.  Z.-Y.;  Li,  X.-D.; 
Xia,  Y.  MRS  Bull.  2005,  30,  356-361. 

(57)  Murphy,  C.  J.;  Sau,  T.  K.;  Gole,  A.  M.;  Orendorff,  C.  J.;  Gao,  J.-X.; 
Gou,  L.-F.;  Hunyadi,  S.  E.;  Li,  T.  J.  Phys.  Chem.  B  2005, 109,  13857- 
13870. 
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Figure  7.  Comparison  of  the  geometry  (a)  of  the  {111}  surface  and  (b)  of 
the  {100}  surface  with  numerical  values  for  Au.  The  surface  atomic  layer 
is  represented  by  spheres  and  subjacent  repeat  layers  as  smaller  spheres 
and  crosses.  Dotted  black  lines  indicate  energetically  favorable  orientations 
of  aromatic  rings  such  as  in  Tyr,  Phe,  and  Trp.  In  (a),  the  ring  or  a  sp2- 
hybridized  group  can  occupy  free  fee  “lattice”  sites,  which  contributes  to 
high  binding  strength.  In  (b),  complementarity  to  chain  molecules  is  difficult 
to  achieve  so  that  small,  flexible  solvent  molecules  and  possibly  designed 
molecules  with  a  rectangular  repeat  are  of  higher  binding  strength. 

charged  (— 3e  per  peptide  plus  3  cations),  which  supports 
precipitation  as  a  result  of  high  local  ionic  concentrations. 
Infrared  spectra  of  the  peptides  in  solution  and  on  the  surface 
exhibit  characteristic  changes  in  the  amide  N— H  stretch,  the 
amide  C=0  stretch,  and  C02  symmetric  stretch  frequencies.7 
We  surmise  that  the  proximity  of  binding  residues  to  the  metal 
surface  and  rearrangements  in  the  local  environment  of  hydrogen 
bonds  (Figures  3  and  4)  account  for  these  observations;  however, 
the  assignment  of  the  superimposed  vibration  frequencies  to 
molecular  details  in  the  simulation  remains  as  a  challenge. 

Hybrid  gold— palladium  nanostructures  were  synthesized 
starting  with  Au  nanoparticles  decorated  with  fusion  peptides 
Fig— Na3~ A3  and  A3— Fig— Na3  and  peptide  A3.  Solutions  of 
the  decorated  Au  nanoparticles  were  exposed  to  a  macroscopic 
Pd  {111}  substrate,8  and  the  affinity  to  the  Pd  {111}  surface 
decreased  in  the  same  order.  The  trend  of  stronger  adsorption 
of  the  peptides  to  Pd  {111}  surfaces  compared  to  that  observed 
on  Au  {111}  surfaces  agrees  with  the  simulation  (Figure  7); 
however,  it  is  unknown  from  experiment  whether  the  N-terminal 
or  the  C-terminal  of  the  fusion  peptide  remains  attached  to  the 
Au  {111}  surface.  The  same  trend  in  affinity  of  the  Au 
nanoparticles  coated  with  Fig— Na3— A3,  A3— Fig— Na3,  and  A3 
was  observed  toward  Pd  nanoparticles  (~3  nm  diameter)  upon 
incubation  with  IGPdClr,  solution  and  reduction  by  NaBH4.8  The 
formation  of  extended  bimetallic  Pd— Au  interfaces8  also 
supports  the  attraction  of  certain  peptides  (A3)  to  bimetal 
interfaces  (Table  4),  although  specific  comparisons  are  not 
possible  on  the  basis  of  available  experimental  data. 

Willett  et  al.17  studied  the  adsorption  of  20  peptides,  each 
comprised  of  8—10  identical  residues  (similar  to  GlylO  and 
Pro  10),  with  fluorescent  tags  at  the  N-terminus  on  various 
surfaces,  including  Au  and  Pd.  Unfortunately,  the  nature  of 
the  surfaces  is  unknown,  and  thermodynamically  preferred 
{111}  surfaces  were  likely  destroyed  by  application  of  a  4-min 
oxygen  plasma  etch  in  favor  of  { 100}  surfaces.  After  incubation 
with  peptide  solutions  at  pH  =  7,  the  number  of  adhered 
peptides  per  surface  area  was  determined  from  the  fluorescence 
intensity  after  several  washing  cycles.  Fluorescence  signals  may 
have  been  quenched  in  peptides  with  all-aromatic  side  chains 
(TyrlO,  PhelO,  etc.),  so  the  interpretation  of  the  data  is  quite 
speculative.  Nevertheless,  strong  adsorption  of  some  peptides 
on  oxidized  Al  and  silica  surfaces  was  found,  while  adsorption 


on  Au  and  Pd  surfaces  was  always  less  than  10%  on  a  relative 
scale:  on  Au,  Arg  and  Thr  were  shown  to  adsorb  slightly,  and 
on  Pd,  His,  Thr,  Tyr,  and  Trp  adsorbed.  Data  from  other 
groups14'7'8  and  simulation  results  leave  no  doubt  that  peptides 
strongly  adsorb  on  Au  {111}  and  Pd  {111}  surfaces,  comparable 
to  their  adsorption  on  surfaces  such  as  alumina  and  silica.  The 
thoroughly  weak  adsorption  of  all  peptides  in  Willett’s  study, 
therefore,  suggests  the  presence  of  {100}  surfaces  on  which 
the  simulation  only  shows  approximately  10%  of  the  adsorption 
energy  compared  to  that  observed  on  {111}  surfaces  (Table  2). 
The  approach  of  Tyr  and  Ser  (similar  to  Thr)  in  A3  to  Au  { 100} 
and  Pd  {100}  surfaces  was  found  in  the  simulation  (Table  2). 

5.2.  Binding  Mechanism.  Overall,  computed  adsorption  ener¬ 
gies  provide  interesting  trends  (Tables  2  and  4)  and  are  in  a 
justified  relationship  with  available  experimental  data.  All 
peptides  containing  8—12  amino  acids,  including  the  nonspe- 
cialized  reference  peptides  GlylO  and  Pro  10,  exhibit  neutral  or 
attractive  interactions  with  gold  and  palladium  surfaces  in  a 
gross  range  from  0  to  —100  keal/mol  of  peptide.  These 
interactions  are  comparable  to  the  strength  of  0—30  hydrogen 
bonds  and  appear  to  be  a  reasonable  order  of  magnitude  for 
metal  binding:  stronger  than  van  der  Waals  interactions  and 
hydrogen  bonds  alone,  and  weaker  than  covalent  interactions.1 
All  residues  or  a  subset  of  binding  residues  in  each  peptide  are 
then  in  direct  contact  with  the  metal  surface  and  account  for 
most  of  the  adsorption  energy  (Tables  2  and  4).  The  amount  of 
attraction  depends  on  the  metal,  the  surface,  the  amino  acid, 
changes  in  conformation,  polarity  of  the  surface  (Pd— Au), 
polarization  at  the  interface,  covalent  contributions  to  bonding, 
and  other  factors.  Direct  contact  with  the  metal  surface  may 
thus  not  be  enough  for  strong  adsorption,  such  as  for  Fig— Na3 
on  Au  {100}  and  Pd  {100}  surfaces. 

The  marked  differences  between  adsorption  on  {111}  and 
{100}  surfaces  can  be  related  to  the  geometry  of  the  surface 
(Figure  7),  in  agreement  with  the  energy  model  of  the  simulation 
(eq  SI,  Supporting  Information).  Adsorbed  molecules  are  most 
strongly  bound  to  the  metal  surface  when  they  serve  as  an 
approximate  extension  of  the  fee  lattice.  This  is  particularly 
evident  on  the  {111}  surface,  where  metal  atoms  in  the  surface 
plane  arrange  in  a  hexagonal  lattice  (Figure  7a).  The  lattice 
involves  angles  of  120°  and  projected  interatomic  distances  of 
1 .66  A  for  Au  and  1 .59  A  for  Pd  in  the  two  atomic  layers  below 
the  surface  atomic  plane  which  create  “free”  fee  lattice  positions 
above  the  metal  surface  with  maximum  van  der  Waals  (metal¬ 
bonding)  interactions.  Therefore,  adsorbed  aromatic  rings  and 
extended  sp2-hybridized  structures  with  bond  angles  of  120° 
and  bond  lengths  of  ~  1 .4 1  A  fit  reasonably  well  on  top  of  the 
surface.  The  density  of  surface  atoms  per  unit  area  is  also 
slightly  higher  on  the  {111}  surface  than  on  the  {100}  surface. 
On  the  other  side,  the  {100}  surface  with  angles  of  90°  and  a 
projected  spacing  between  “empty”  fee  lattice  sites  of  2.88  A 
for  Au  and  of  2.75  A  for  Pd  is  much  less  complementary  to 
common  sp2-  or  sp3-hybridized  peptide  chains  or  other  organic 
molecules  (Figure  7b).  This  quadratic  pattern  of  available  fee 
lattice  sites  induces  a  significant  mismatch  to  usual  bond  lengths 
and  bond  angles,  which  makes  it  difficult  to  bind  peptides  or 
other  chain  molecules  to  { 100}  surfaces.  Instead,  small,  mobile 
water  molecules  are  preferentially  adsorbed  and  lead  to  small 
adsorption  energies  for  all  peptides  on  {100}  surfaces  (Table 
2).  In  summary,  aromatic  residues,  sp2-configured  groups  such 
as  C02  of  Asp,  and  sp3-hybridized  alkyl  chains  with  tetrahedral 
bond  angles  of  109°  and  bond  lengths  of  1.54  A  are  thus  a  better 
fit  to  {111}  surfaces  compared  to  {100}  surfaces  (Figures  3,  4, 
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and  6),  and  the  compatibility  with  angles  and  spacing  between 
available  lattice  sites  can  be  exploited  to  design  strongly  or 
weakly  binding  peptides. 

The  nature  of  noncovalent  binding  of  peptides  and  organic 
molecules  to  metal  surfaces  may  thus  be  described  as  a  soft 
epitaxy  driven  by  complementarity  to  “empty”  lattice  positions. 
Yet  bound  molecules  possess  significant  translational  freedom 
compared  to  atoms  in  the  metal  surface  plane,  so  several  factors 
contribute  to  binding  strength:  (1)  the  interaction  strength  of 
the  metal,  proportional  to  its  surface  energy;  (2)  the  interaction 
strength  of  the  peptide  or  individual  residues  and  its  comple¬ 
mentarity  to  the  metal  surface;  (3)  the  interaction  strength  of 
the  solvent  and  its  complementarity  to  the  metal  surface;  and 
(4)  entropic  factors  and  changes  in  secondary  structure  upon 
adsorption.  In  addition,  polarization  of  the  metal  surface  adds 
an  amount  of  attractive  energy  that  would  be  similar  on  {111}, 
{ 100},  { 1 10},  or  other  surfaces.  Also,  charge  transfer  is  possible. 
Both  of  these  effects  will  be  less  significant  in  the  presence  of 
a  solvent  (e.g.,  water)  of  polarity  similar  to  or  higher  than  that 
of  the  peptide;  however,  they  may  tune  peptide  adsorption  on 
less  attractive  (e.g.,  {100})  surfaces. 

The  concept  of  soft  epitaxy  also  suggests  that  peptide 
adsorption  is  not  highly  specific  to  individual  amino  acids  but 
rather  to  the  peptide  as  a  whole,  and  this  aids  in  understanding 
the  self-assembly  of  larger  structures,  such  as  Met  nanogratings 
on  even  Ag  {111}  surfaces.12  Additional  evidence  for  the 
proposed  binding  mechanism  from  earlier  simulations  is  pro¬ 
vided  in  section  S3  of  the  Supporting  Information. 

It  is  also  instructive  to  compare  the  noncovalent  mechanism 
of  peptide  adsorption  to  the  covalent-type  adsorption  of  al- 
kanethiols  on  Au  surfaces. 14-16,58  The  dissociation  energy  of 
strongly  binding  peptide  groups  such  as  C02  (Asp)  and  PhOH 
(Tyr)  is  on  the  order  of  10—20  kcal/mol  (sections  3  and  4), 
which  is  similar  to  an  estimated  dissociation  energy  of  11.5 
kcal/mol58  of  the  Au— S(C„H2jj+i)  bond  in  thiols,  although  the 
value  for  thiols  still  remains  uncertain.  A  major  difference 
between  peptides  and  thiols  is  that  peptides  (excluding  Cys) 
spread  out  their  interaction  over  a  greater  area  and  several 
functional  groups,  while  alkanethiols  lead  to  pinpoint  interac¬ 
tions  at  the  R— SH  headgroup  which  induce  locally  strong 
covalent  forces.  Thus,  alkanethiols  are  likely  to  support  the 
formation  of  adatoms  and  vacancies  during  the  reduction  of 
precursor  metal  salts.1416  Besides,  the  high  density  of  Au  atoms 
on  even  surfaces  (three  Au  atoms  per  thiol  headgroup)15  favors 
thiol  binding  when  the  local  surface  roughness  (i.e.,  the  effective 
surface  area)  is  increased,  while  the  soft  epitaxial  adsorption  of 
peptides  benefits  from  locally  regular,  complementary  surfaces. 

6.  Conclusion 

Using  molecular  dynamics  simulation  and  the  consistent 
valence  force  field  with  accurate  extensions  for  fee  metals,  we 
computed  adsorption  energies  of  strongly  binding  peptides 
identified  by  phage  display  (A3,  Fig— Na3,  Pd2— Cl,  Pd4— Cl) 
as  well  as  control  peptides  (Flgd— Na2,  GlylO,  ProlO)  on  {111} 
and  {100}  surfaces  of  Au  and  Pd  and  on  Pd— Au  {111}  bimetal 
surfaces  in  aqueous  solution.  Adsorption  energies  range  from 
0  to  —100  kcal/mol  per  peptide  of  8—12  amino  acids.  Strongly 
binding  peptides  often  assume  a  flat-on  2D  conformation  so 
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that  most  residues  are  in  direct  contact  with  the  metal  surface, 
and  weakly  binding  peptides  retain  more  features  of  their 
preferred  3D  conformation  in  solution  so  that  layers  of  water 
molecules  are  present  between  the  surface  and  a  number  of 
peptide  residues. 

A  clear  preference  for  adsorption  to  {111}  surfaces  over 
{100}  surfaces  is  seen  for  all  peptides  in  aqueous  solution,  which 
can  be  explained  under  consideration  of  available  fee  lattice 
sites  above  the  metal  substrate.  Above  the  {111}  surface  of  Au 
and  Pd,  a  hexagonal  array  of  available  lattice  sites  with  a  spacing 
near  1.60  A  can  be  favorably  populated  by  aromatic  rings  or 
sp2-hybridized  groups  in  the  side  chain  or  by  the  main  chain  of 
the  peptide.  On  the  other  hand,  a  rectangular  pattern  of  available 
lattice  sites  with  a  spacing  of  2.88  A  above  the  {100}  surface 
is  more  difficult  to  populate  on  time  average  by  peptide 
molecules  so  that  direct  contact  with  small  mobile  water 
molecules  is  preferred.  Besides,  all  peptides  under  consideration 
indicate  a  stronger  adsorption  on  Pd  {111}  surfaces  compared 
to  Au  {111}  surfaces,  which  appears  to  be  related  to  the  higher 
surface  tension  and  thus  higher  attraction  to  Pd  {111}  surfaces 
compared  to  Au  {111}  surfaces.  Pd— Au  {111}  bimetallic 
surfaces  also  possess  additional  polarity,  approximated  by  atomic 
charges  of  +0.3e  and  — 0.3e  at  the  Pd  and  Au  sides  of  the 
interface,  which  influences  the  orientation  of  polar  groups  in 
the  interfacial  region,  such  as  —  C(=0)— NH—  peptide  bonds, 
NH3+  groups  in  Lys,  or  C02-  groups  in  Asp.  Adsorption 
energies  are  found  to  be  on  the  order  of  10  kcal/mol  higher 
(A3)  and  lower  (Fig— Na3)  compared  to  those  observed  on  the 
most  attractive  monometallic  surface. 

The  concept  of  complementarity  to  available  fee  lattice 
positions  allows  the  distinction  of  strongly  binding  residues  and 
less-binding  residues.  Strong  adsorption  on  Au  {111}  and  Pd 
{111}  surfaces  is  particularly  supported  by  aromatic  and 
conjugated  residues  such  as  Tyr,  Phe  (preferably  on  Pd),  His, 
Arg,  Asp,  and  Asn,  which  can  be  regarded  as  “specific  binding”. 
However,  strong  adsoiption  also  seems  possible  by  less  attractive 
residues  such  as  Gly,  Thr,  Ser,  and  others.  The  binding  is  then 
enhanced  through  cooperative  effects  between  many  residues 
along  the  chain,  e.g.,  GlylO,  and  may  be  regarded  as  “nonspe¬ 
cific  binding”. 

The  simulation  approach  presented  is  a  first  attempt  to  model 
metal  interfaces  with  water  and  (bio)organic  moieties  using 
validated  force  fields,24  explicit  water  models,  and  physically 
justifiable  peptide  concentrations  (~50  mM)  at  pH  =  7  on  a 
time  scale  of  10  ns.  We  employed  a  simple  and  efficient  method 
to  compute  energies,  free  energies,  and  entropies  of  adsorption 
of  solute  molecules  on  a  given  surface  that  involves  the 
simulation  of  four  systems,  the  surface— solute— solvent  system, 
the  solute— solvent  system,  the  solvent  system,  and  the 
surface— solvent  system,  under  consideration  of  equal  molecular 
volumes  of  each  component  under  NVT  conditions.  Limitations 
of  the  semiempirical  classical  simulation  include  the  neglect  of 
the  explicit  electronic  structure,  approximations  in  the  energy 
model  as  well  as  a  limited  accuracy  of  the  force  field 
(particularly  for  water  and  peptides),  and  the  neglect  of 
polarization  effects  (of  limited  importance  on  even  surfaces). 
Boltzmann  averages  of  the  sometimes  strongly  adsorbed  pep¬ 
tides  are  challenging  to  obtain  in  spite  of  several  independent 
start  conformations,  and  the  finite  size  of  3—5  nm  in  the  3D 
periodic  model  imposes  a  small  but  artificial  symmetry  restraint. 

The  results  are  intended  as  a  first  guide  toward  understanding 
the  binding  affinities  of  peptides  to  metal  and  bimetal  surfaces, 
aiming  at  the  design  of  peptides  for  specific  surfaces.  A  range 
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of  experimental  data  supports  the  concept  of  dynamic  epitaxy 
in  solution  and  the  description  of  binding  versus  less-binding 
residues,  as  well  as  computed  binding  energies  in  semiquanti- 
tative  accuracy  or  better.  Also,  quantum-mechanical  studies  on 
peptide  fragments  support  our  view.  Nevertheless,  significant 
further  insight  could  be  obtained  through  further  specific 
experimental  characterization  of  binding  constants  and  NMR 
data  and  through  the  exploration  of  the  nanomechanical 
properties.55'56  Other  future  challenges  include  the  quantitative 
evaluation  of  the  effects  of  curved  and  stepped  surfaces,  irregular 
surface  structures,  and  the  presence  of  ions  on  the  binding  of 
peptides,  as  well  as  testing  of  specifically  synthesized  molecules 
for  complementarity  to  {111}  and  {100}  surfaces. 
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